Magnon-mediated interactions between fermions depend strongly on the lattice 

structure 
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We propose two new methods to calculate exactly the spectrum of two spin-i charge carriers 
moving in a ferromagnetic background, at zero temperature. We find that if the spins are located on 
a different sublattice than that on which the fermions move, magnon-mediated effective interactions 
are very strong and can bind the fermions into low-energy bipolarons with triplet character. This 
never happens in models where spins and charge carriers share the same lattice, whether they are 
in the same band or in different bands. This proves that effective one-lattice models do not describe 
correctly the low-energy part of the two-carrier spectrum of a two-sublattice model, even though 
they may describe the low-energy single-carrier spectrum appropriately. 

PACS numbers: 71.10.Fd, 71.27. +a, 75.50.Dd 



When studying properties of complex materials with 
ions on several sub-lattices it is customary to use sim- 
plified, one-lattice Hamiltonians to describe their low- 
energy physics. For example, instead of a two-sublattice 
model for the CuC>2 plane including both Cu 3d x 2_ y 2 and 
O 2p orbitals [1], many prefer a one-band Hubbard model 
on the Cu (sub)lattice. States in this simpler model are 
Zhang- Rice singlets (ZRS), i.e. bound singlets between 
a hole at a Cu site and a doping hole delocalized over its 
four O neighbours with the same d x i_ y i symmetry [2]. 

Such a composite object may describe well the low- 
energy quasiparticle, although this is still debated [3]. 
Less clear is whether a model based on such states that 
mix together charge and spin degrees of freedom, can 
properly describe quasiparticle interactions, especially 
those mediated through spin fluctuations. Most oxides 
have at least one phase with long-range magnetic order, 
and magnon exchange is believed by some to be a key 
component determining their properties, eg. as the main 
"glue" for pairing in cuprates, which likely controls the 
value of T c [4]. 

Here we show that effective one-lattice models severely 
underestimate the magnon-mediated attraction between 
carriers, compared to their two-sublattice "parent" 
model. The magnetic background is chosen as ferro- 
magnetic (FM). This is much simpler than an antifer- 
romagnetic (AFM) background, but it allows for exact 
solutions. Thus, any qualitative differences are inher- 
ent to the models themselves. Moreover, our conclusions 
are relevant to the modeling of carriers in AFM back- 
grounds, and raise serious questions about the ability of 
ZRS-likc constructs to describe correctly low-energy two- 
carrier states. 

Our models are sketched in Fig. 1. For simplicity, 
we address the one-dimensional (ID) case; generaliza- 
tions are straightforward. Model I is the "parent" , two- 
sublattice, two-band model. Model II is a two-band, sin- 
gle lattice effective model. Model III is an even simpler 
one-band effective model. 



In models I and II, one band hosts the spin S degrees 
of freedom, described by: 



n s = -j 



Si ■ Si+i — s 



This favours a FM ground-state \FM) = \ + S, . . . , +S) 
for the undoped system. Spin-i doping charge carriers 
occupy states in another band. In model I, this is located 
on a different sublattice, for example like in a CuO chain 
with spins on Cu and holes on O sites. In model II, 
they are on the same lattice. In both cases carriers are 
described by a Hubbard model: 



He = -tJ2 (4+l+5,a C i+6,° + hx ) + U J2 fl 



where 5 = \ for model I, 5 = for model II, Ci+s,a is 
the annihilation operator for a carrier of spin a from site 
% + 5, and n i+s ^ = c\ +S a c i+s ^. 

Interactions between charges and spins are described 
by the simplest exchange model. In model I, a carrier 
interacts with its two neighbour spins: 
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FIG. 1. Models I and II have two bands: one occupied by 
spins (arrows), and one (empty circles) hosting carriers intro- 
duced by doping (filled circles, with arrow showing the spin). 
In the "parent" model I, these are on different sublattices. In 
model II, they are on the same lattice. Model III has one 
band which hosts both spins (arrows) and ZRS-like polaron 
cores (filled circle). 
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while for model II there is on-site exchange: 

i 

Here s i+s = \ Y, a .j3 4+S, a VapCi+8,p, where a are Pauli 
matrices. We set h = 1 and the lattice constant a = 1. 

The one-band model III is described by the Hubbard 
Hamiltonian Ti c projected onto the appropriate subspace, 
eg. St ot = NS — h if one spin-down carrier is added. 

We set S = | (higher S is discussed elsewhere [5]), 
and first review the one-carrier case. If a =t, the T = 
problem is trivial in models I and II: spin-flips are im- 
possible, therefore the eigenstates c\AFM) have energy 
Ek\ = £k + jJo, where 7 = 5 (7 = j) in model I (II), and 
€k = — 2tcos(k) is the free carrier dispersion. Model III 
has a serious problem: it cannot distinguish an undoped 
system from one doped with spin-up carriers. 

The interesting case is for a carrier injected with spin- 
down. In model III, this is trivial: the "hole" moves freely 
with energy ek ■ The exact solution (in any dimension) for 
the Green's function G i (fc,w) = (FM \c k> x.G(w)c\^\F M) 

where G(u>) = [u> + irj — H}^ 1 is the resolvent for H = 
H s +H c + 'H ex , has long been known for model II [6]. Its 
generalization to two-sublattice models was proposed re- 
cently [7]. The off-diagonal part of H ex mixes c^JFM), 
of energy = e k — 7 Jo, with the continuum of one- 
magnon states c\_ q1 .S q '\FM) of energy E^-g^ + fi g . 
Here, S~ — -7= J2i e lqRi S^ is the magnon creation op- 
erator and fig = 2Jsin 2 (g/2) is the spin- wave dispersion. 
Also, S~ = Sf — iSf is the spin-lowering operator and 
N — > 00 is the number of sites on either sublatticc. 

The interesting case is Jo > 0: hybridization pushes 
the discrete state further below the continuum, result- 
ing in a low-energy infinitely lived quasiparticle (spin- 
polaron), discussed next. If J < 0, the low energy states 
are the incoherent continuum describing the scattering 
of the carrier off a free magnon [6] . One cannot further 
simplify the description of such states. 

Returning to the spin-polaron that is the low-energy 
quasiparticle for J > 0, its structure can be understood 
in the limit Jo ^ t, J- We start with model II. T-lix 1 ^ is 
minimized by an on-site singlet between the carrier and 
its spin, \s)i = ^cjj_ — c\^S~^j \FM), with all other 
spins in the FM state. The energy of this degenerate 
state is Hi I x I ^\s)i = — |j |s};. Hopping lifts degeneracy, 
and to first order in t, J, the polaron energy is: 

E^\k)^-\.h + \e k+ J - (1) 

with \P\i,k) = -^^2i etkRi \ s )i- Thus, the spin-polaron 
is an on-site singlet between the charge and its local spin 
(or a bound state of the carrier with a magnon at the 
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FIG. 2. (color online) (a) Model I, and (b) model II density of 
states pi(k,cj) — — ^ImG^A:, oj). Contour plots show exact 
results. Full lines are Eqs. (1), (2), while dashed lines mark 
the expected onset of the continuum, at min, (Ek- q ,f + fig). 
Here J/t = 0.05, J /t = 5, <q/t = 0.02. 

same site) that moves with an effective hopping t/2 sup- 
pressed by the magnon cloud overlap. The last term is 
the FM exchange energy lost in the magnon's presence. 

Similar considerations apply to model I. Again, we 
only discuss the case Jo > which has a low-energy 
quasiparticle. Because of the two-sublattice structure, 
the ground state of Hix is the three-spin polaron (3sp) 

|3 SP ) i+ , = (v^W^U/r^ 1 ) ' FM >< of en " 

ergy %i^|3sp) l+ i = -J |3.sp) l+ i [8]. It describes a 
bound state between the carrier and a magnon on either 
neighbouring spin. Hopping lifts degeneracy, resulting in 
\Pi, k ) = VW e^ R *+^\3sp) l+ i with energy: 

4 /) (fc)«-J + jU + ^. (2) 

In Fig. 2 we compare the exact spectra, calculated as 
in Refs. [6, 7], with these asymptotic expressions (lines) 
for models I and II. The agreement is excellent for large 
Jo and remains reasonable even down to Jo ~ t, showing 
that this singlet/3sp description of the spin-polaron is 
robust. Fig. 2(a) also shows another quasiparticle below 
the continuum. This is based on a higher eigenstate of 
Hex (for more details, see Ref. [7]. Here we focus only 
on the low-energy physics) . 

Based on this, it is clear that if we are only 
interested in the low-energy spin-down quasiparti- 
cle, we can map model I onto II and III, after 
some rescaling of parameters. The mapping from 
model I to II is obvious if we rewrite \Pi,k) — 
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a ikR d J_ 
V2 



dl 



l-FM), i.e. a singlet- 



's cl 



type state like |Pn,fc), except that the carrier is in the 
"on-site" orbital dl , = 4= [e'^ct, i H 

Thus, |Pi,fc) is the direct analog of the ZRS: in the 
k = ground-state, the carrier occupies a linear combi- 
nation of neighbour orbitals with the same s-type sym- 
metry as the orbital giving rise to the Si spin, with which 
it locks in a singlet. (If the orbitals had other symme- 
tries, the hopping matrices would be different and the 
phases would change accordingly [7]). Like in the ZRS, 
d\ irT is not orthogonal to its neighbour d\ i±1(T . Unlike 
the ZRS, these "on-site" states depend on k. This en- 
sures the normalization of \Pi,k) in the entire Brillouin 
zone, unlike for the ZRS whose normalization diverges at 
k = [2]. 

Nevertheless, if we ignore such complications and re- 
place dfc,i,<7 - >■ Cio-, model I maps onto model II as far 
as the low-energy quasiparticle is concerned. Mapping 
to model III is the next step of replacing the ZRS with 
a "hole" which lives in the same band as the spins of 
the undoped system. Since model III has a quasiparticle 
band in the S^ ot = (N — l)g sector, it can also be mapped 
onto the quasiparticle band of model I, after rescaling. 

The important question is whether this mapping be- 
tween low-energy sectors carries on to cases with more 
carriers. We consider the two-carriers case, and show 
that model I has low-energy states which have no equiv- 
alent in models II and III. Specifically, magnon-mediated 
interactions can stabilize a low-energy bound pair in 
model I, if Jo is sufficiently large. Such states never occur 
at low energies in model II, and are impossible in model 
III. This proves that modeling the proper lattice struc- 
ture is essential to correctly describe magnon-mediated 
interactions. 

The two-particle problem also has several cases. The 
trivial one is for two spin-up carriers. In models I and 
II the carriers behave like non-interacting particles, with 
eigenstates c^cJ./J-FM) of energy Ek^+Ek',f Model III 
does not distinguish this case from an undoped system. 

The most interesting case is when one carrier is injected 
with spin up and the other with spin down. We analyze 
it in detail for models I and II, and then briefly discuss 
the case where both carriers are injected with spin-down. 
We solve the problem exactly using two new methods. 

The first method calculates two-particle Green's func- 
tions in momentum space. Because of translational in- 
variance, total momentum is conserved so non- vanishing 
matrix elements are G(k,q,q' ,uj) = (k,q\G(ui)\k,q'), 
where \k, q) = c\ + ^ ct _ q JFM) is a two-particle state 

with total momentum k. Momentum q is not conserved: 
on-site scattering changes it, as do magnon-mediated in- 
teractions. In the latter, the spin-down carrier flips its 
spin creating a magnon, followed by absorption of the 
magnon by the other (initially spin-up) carrier. The 




0.4 



0.3 



0.2 



< 0.1 
0.0 





-10 


co/t 


10 


20 


1 


1 


II 




i 




continuum 


(b) 
xiot 




continuum 


- 




y — \ 
/ U=0\ 


s 




> 




xlO 2 


ft: 




- 








/ i , V 




-10 


0, 

co/t 




10 



FIG. 3. (color online) Spectral weight A(k = 0, n = n' = 1, w) 
for model I (a) and II (b). Expected continua locations are 
marked, as are triplet (T) and singlet (S) bipolarons (arrows). 
Here J /t = 20, J/t = 0.05, rj/t = 0.1 and U/t = 0, 1,5. 



magnon's momentum is, thus, transferred from one car- 
rier to the other, and their spins are exchanged. It is 
precisely this effective interaction that interests us. 

Because only one-magnon states are accessible from 
the original state, we obtain a closed system of equations 
of motion for these Green's functions. This is rather sim- 
ilar to the single-carrier case [6, 7], however the solution 
is now less trivial. In the Supplemental Material [9] we 
present the steps to reduce it to a closed equation that 
can be efficiently solved for a finite chain with N ~ 100. 
However, this becomes costly in higher dimensions, and 
may not generalize to other cases (eg., more carriers). 

We also propose a real-space solution for the infinite 
chain which allows us to find the symmetry of the 
pair (singlet vs. triplet), and also generalizes to more 
carriers. It is based on the few-particle solution of Rcf. 
[10]. Our case is more complex because when present, 
the magnon also counts as a "particle" . Thus, the 
system switches between having the carriers in states 



like 



and three "particle" configurations \k,n,m) = 
_^ e^+^)c\ +s ^ l+s+n ^S- +m \FM). Exchange 

connects propagators G(k, n, n' , uS) = (k, n\G(cj)\k, ml) to 
G(k,n,n' ,m,uj) = (k,n\G(ui)\k,n',m), and vice-versa. 
The solution of these equations is described in [9] . 

Poles of these propagators mark the two-carrier spec- 
trum. Fig. 3 shows A(k, n, n', u>) = — ImG(fc,n,n',w) 
for k = 0,n = n' = 1 and Jo = 20t, chosen so large to 
simplify the task of identifying features. 

The expected features in the spectrum of model I are: 
(i) a continuum describing states where the spin-down 
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FIG. 4. (color online) The three highest weight configurations 
contributing to the low-energy bipolaron of model I. 



carrier forms a polaron of momentum k! and the spin- 
up carrier has momentum k — k' . These states span 
{Ep\k') + Ek-k',-f}k' so the band-edges are known from 
single-carrier results. Using Eq. (2) and ignoring small J 
corrections, the band edges are here expected at — 13. 7 1 
and — 6.3t (dashed lines), in good agreement with the re- 
sults; (ii) a continuum {E^i^ + E^-^i-q^ + de- 
scribing states comprised of two spin-up carriers and 
a free magnon. Expected band-edges at l&t and 24i 
(dashed lines) show again good agreement; (iii) if there 
is a higher energy polaron state, like in Fig. 2(a), it will 
also generate a continuum. Here, its edges should be at 
9t and lit. It is not seen in Fig. 3(a) because of vanish- 
ing overlap with the \k = 0,n= 1) state, but it appears 
in momentum-space Green's functions [9]. 

Any states outside these continua are bound states, 
"glued" through magnon-exchange. Five such bipolarons 
appear in Fig. 3(a) (the apparent overlap between SI 
and S3 with their nearby continua is due to the broad- 
ening rj = O.lt). For k = 0, the real-space solution 
allows us to identify two as triplets and the other as 
singlets (see arrows). For finite k, there is no def- 
inite symmetry: these bound states have finite over- 
lap with both singlet and triplet configurations, such as 



Consider now the low-energy bipolaron. Its largest 
overlap is with the configurations sketched in Fig. 4, 
hybridized through the Jo exchange. It is the ability 
of both carriers to interact with the same spin to ex- 
change the magnon, that stabilizes this state. This also 
explains why it is a triplet at k = (symmetric com- 
bination is favoured), and its insensitivity to U: unlike 
singlets, triplets do not permit double occupancy. This 
bipolaron forms if J /t > 6.5, and is quite light [5]. 

Model II has different behaviour, see Fig. 3(b). Again, 
a polaron+spin-up carrier continuum, and a magnon + 
two spin-up carriers continuum, are expected and ob- 
served. If Jo is sufficiently large that the two do not 
overlap, a bound state appears between them, at uj ~ U. 
This bipolaron is a singlet at k = 0, quite similar to 
the S2 state in model I. It has very little weight in the 
\k = 0,n = 1) configuration, which is why we had to 
magnify the peaks. Most of its weight is on the on-site 
(n = 0) configuration [9]. No triplet-like bound states 
appear, and no low-energy bipolaron is possible for any 
values of the parameters. This is not so surprising, since 
in this model only one carrier can interact with a given 
spin/magnon. If the carriers are on the same site they 



form a singlet which has no interactions with the local 
spin. If they are on neighbouring sites, magnon exchange 
(now controlled by J, not by Jo as in Fig. 4) is too weak 
to stabilize a low-energy bipolaron. Thus, model II sim- 
ply cannot describe the low-energy physics of model I 
in this sector. Model III fails as well, since it does not 
distinguish between a spin-up carrier and a lattice spin. 

We also considered the two spin-down carrier states. 
For models I and II, the t = solutions suggest that the 
lowest energy feature is the two-polaron continuum. This 
is reasonable, since each carrier can bind its own magnon 
to create a polaron. Simultaneous interaction of one car- 
rier with two magnons is impossible in model II, and 
while possible, it is energetically unfavourable in model 
I, so low-energy bipolarons do not appear in this sector. 
Turning on hopping further favours the continuum, since 
polarons are lighter than a bipolaron [5]. In model III, 
no interactions are possible between two "holes" , since 
they simply reshuffle the FM spins as they move. This is 
why we have not investigated this case in more detail. 

To summarize, we extended the exact solution for a 
single charge in a FM background to cases with two or 
more carriers, and discussed in detail the nontrivial case 
where one carrier is injected with spin up and the other 
with spin down. The low-energy physics depends essen- 
tially on the model. If the spins are intercalated between 
the carrier sites, magnon-exchange is enhanced and can 
bind low-energy bipolarons. If spins and carriers live on 
the same lattice, such low-energy states are impossible. 
If spins and carriers live in the same band, there are 
no magnon-exchange interactions in any of the allowed 
cases. These three models have very different low-energy 
states in the two-particle sector, even though their one- 
particle sectors can be mapped onto one another. 

This shows that in order to properly describe magnon- 
mediated interactions, one must use the proper sublatticc 
structure, not simpler effective one-lattice models. While 
our work is for a FM background, it is directly relevant 
for AFM backgrounds as well, since the magnon exchange 
is a rather local process. Based on our results, for a Cu02 
lattice one should expect strong magnon-mediated inter- 
actions between two holes located at neighbour O sites, 
through their common Cu (the 2D analog of Fig. 4). A 
one-band model based on ZRSs simply cannot describe 
this process, and is therefore likely to severely underesti- 
mate the role of magnons as the "glue" for pairing. 
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